FONCTIONS DE MITTAG-LEFFLER ET PROCESSUS DE LÉVY STABLES 

SANS SAUT NÉGATIF 



THOMAS SIMON 



RÉSUMÉ. On remarque que la fonction E a (x a ) — ax a ~ x E' a {x a ) est complètement monotone pour 
tout a € [1,2]. Grâce à l'expression de sa densité de Bernstein, on en déduit une identité en loi entre 
suprema de processus de Lévy a— stables complètement asymétriques. Dans le cas spectralement 
positif, on retrouve l'expression d'une constante de petites déviations unilatères qui avait été obtenue 
avec une autre méthode par Bernyk, Dalang et Peskir pQ. 



1. Une propriété des fonctions de Mittag-Leffler 

La fonction de Mittag Leffler E a d'indice a > est définie pour tout z G C par le développement 
en série entière 



(i) W = Efn 



,..-,- + an) 

où r désigne la fonction Gamma d'Euler. En écrivant 1/T comme une intégrale curviligne le long 
d'un chemin de Hankel - voir par exemple [3] Chapitre IX. 4 et Figure 56 : 



1 1 



r an - l édt, 



T(l + an) 2vri J H 
et en intégrant la série terme à terme on trouve la représentation bien connue 

(2) EJz) = — / -dt 

K ' K ' 2vri J Hz t a -z 

où H z désigne un chemin de Hankel contournant le cercle centré en l'origine de rayon \z\. En 
particulier on a pour tout x > 

EJx a ) = — / —dt. 



On peut appliquer la même procédure à la fonction suivante définie sur (0, +oo) : 

00 _em-l 
a—1 771/ / a\ \ r x 

^ r an 

n=l 

pour trouver la représentation 

ax a - l E'Jx a ) = — / -dt. 

aK ' 27T1 J„ a t<* - X« 
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On définit maintenant la fonction 



1 f , t a ~ L -x 



KX—1 ~a— 1 



DJx) = EJx a ) - ax^E'Jx 01 ) = / e* dt 

pour tout x > et on la prolonge par continuité en zéro en posant D a (0) = 1. On rappelle qu'une 
fonction / : (0, +oo) — * M. est complètement monotone - ce que nous noterons par CM - si elle est 
indéfiniment dérivable et si 

(-i) n f^(x) > o 

pour tout n G N, x > 0. Le théorème de Bernstein - voir par exemple [5] Chapitre XIII. 4 - dit que 
si / est CM et si f(x)—*l quand x — ► 0, alors il existe une unique mesure de probabilité p sur M + 
telle que 

/■OO 

f(x) = / e- xt n(dt) 



pour tout x > 0, la réciproque étant immédiate. 
Théorème 1. Pour tout a G [1,2], la fonction D a est CM. 

Preuve : La propriété est évidente pour a = 1 puisque D\(x) = et pour a = 2 puisque l'on 
voit facilement que D2{x) = e~ x , prototype d'une fonction CM avec mesure de Bernstein fj, = 5\. 
Supposons maintenant a G (1, 2) et fixons x > 0. En raisonnant comme dans [1] Chapitre VIL 7, la 
fonction 

t a-l _ x a-l 

q : 1 1— > 

i Q - £ a 

se prolonge analytiquement dans tout le plan fendu C/(— oo,0] en ayant posé g(x a ) = (1 — l/a)/x. 
Par le théorème de Cauchy, on a donc 

où H est un chemin de Hankel entourant l'origine indépendamment de x. Pour tout p > 0, on 
décompose classiquement H en la branche inférieure de M~ parcourue de — oo à — p, le cercle centré 
en l'origine de rayon p parcouru dans le sens trigonométrique, et la branche supérieure de M~ 
parcourue de — p à — oo. Comme x > 0, l'intégrale le long du cercle tend vers quand p — > 0. 
L'intégrale le long des deux branches vaut d'autre part 

-in(a— l) s a— 1 _ x a— 1 \ roo I „m(oc— 1) s a— 1 _ x a— 1 \ \ 

; \ ds - I e~ s : ds 




p 

giTra^o— 1 _|_ 1 g— motga— 1 _|_ j-Q— 1 



ds 



gi7ro^a x a g—mctgOt x ct 

quand p — > 0. Après changement de variable s = ra, on trouve 
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1 r°° fe ma u a ~ l + 1 e-' ma u a ~ l + 1 



DJx) = / e~ xu —. ■ : — du 

y ' 2vri ./o V e ma u a - 1 e~ ma u a - 1 ' 



sin7ra f°° _ xu ( u a i (l+u) 



7T Jq \u 2a — 2u a cos ira + 1 

ce qui entraîne que D a est CM avec pour mesure de Bernstein 

3 ^ (dt) = \ +2a ' Z~n 

■n{t za — 2t a cos ira + 1) 



du, 



□ 



Remarques 2. (a) La question de la monotonicité complète ne se pose évidemment pas pour E a 
elle-même mais il est classique [8] que x i— ► E a {—x) est CM pour tout a G (0, 1]. Cette propriété 
avait été obtenue auparavant par Feller avec un argument probabiliste et nous y reviendrons dans 
la section suivante. Nous renvoyons à [TT] pour des résultats plus généraux, qui entraînent en 
particulier que x \—* E a {—x) n'est pas CM quand a > 1. 

(b) On voit facilement que fj, a =3- quand a — » 1 et on peut montrer que fx a ôi quand a — » 2. 
Jean- François Burnol m'a aussi fait remarquer que la fonction a i— ► E a (x) est en fait analytique 
dans un cône contenant IR + pour tout x > 0. Quand a g]0, 1[ les calculs précédents restent valables 
et entraînent que — D a est CM avec mesure de Bernstein 

(sinvrQ)t a - 1 (l + t) 
vr(t 2a - 2i a cosvtq + 1) 

En revanche, quand a > 2 il semble que l'on perde toute monotonicité puisque la fonction g ci-dessus 
admet alors au moins deux pôles à l'intérieur du contour. On peut par exemple calculer 

Di{x) = -(e _a: + cosx + sinx). 

Il n'est pas impossible qu'il faille retrancher ou ajouter des dérivées successives à D a pour retomber 
sur une fonction CM. 

Considérons maintenant pour tout a G [1,2] la fonction F a définie sur M + par 

F a (x) = D^x 1 ^) = E a {x)-ax l - l ' a E' a (x). 

Comme x i— ► x l / a est une fonction positive de dérivée CM, par le Théorème 1 on sait en appliquant 
le Critère 2 dans [5] Chapitre XIII. 4 que F a est également CM. En fait on aurait pu obtenir ce 
résultat, qui est légèrement plus faible que le Théorème 1, par un argument probabiliste. C'est ce 
que nous allons commencer par expliquer dans la section suivante. 
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2. Une identité en loi pour les processus stables 

On fixe a G (1,2] et on considère un processus de Lévy a— stable sans saut négatif X = 
{Xt, t > 0} , normalisé tel que E[e _Xl ] = e. On pose X = —X qui est un processus de Lévy 
a— stable sans saut positif et l'on renvoie aux chapitres VII et VIII de [2] pour plus de détails 
concernant X. La formule de Lévy-Khintchine s'écrit 

E[e- XXt ] = E[e xjit ] = e txa 

pour tout t, A > 0, de sorte qu'avec cette normalisation, quand a = 2 on a X = X = V2B où B 
est un mouvement brownien standard. On définit les suprema 

St = sup{X s , s < t} et St = sup{X s , s < t}, t > 0. 

Pour tout q > 0, soit T q une variable exponentielle de paramètre q indépendante de X. Comme 
St = — inî{X s , s < t} pour tout t, la factorisation de Wiener-Hopf pour X - voir [2] Chapitre VII. 1 
Formule (3) - entraîne que pour tout A > q x l a > 

(4) E[e-^] = Q ^7 ql,a \ 

et après intégration par parties, il vient 

A a-1 1-1A* 



roo 

/ e~ Xx F[S Tq > x] dx 
Jo 



X a — q X a -q 

pour tout À > q l l a > 0. En fait les deux formules ci-dessus sont valables pour tout À, q > avec 
un prolongement par continuité immédiat en À = q l l a . 

Le lien avec les fonctions de Mittag-Leffler est le suivant : en intégrant ([T]) terme à terme on 
peut aussi obtenir l'expression de la transformée de Laplace de E a (qx a ) : pour tout À > q x l a on a 

(5) / e~ Xx E a (qx a ) dx 



X a -q 

et en intégrant par parties, il vient 

X a-1 l-l/a 



- Xx (E a (qx a ) - aiqx^^^E'^qx )) dx 



X a - q X e 

roo 



= / e- Xx ¥[S Tq > x] dx. 
Jo 

Signalons que la formule © avait déjà été utilisée dans pour obtenir de nouvelles correspondances 
symboliques entre fonctions E a d'indices différents, et plus récemment dans [3] pour calculer une 
constante de petites déviations en norme uniforme pour X. L'inversion de la transformée de Laplace 
donne 

(6) nSr q > x] = E a {qx a ) - a(qx a ) 1 - 1 / a E' a (qx a ) 
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pour tout q, x > 0. Il n'est pas évident a priori que la fonction x i— > ^[S T1 > x] soit complètement 
monotone. En revanche, on peut retrouver à l'aide de © le résultat énoncé à la fin de la première 
section. Pour cela, on introduit pour tout x > le premier temps de sortie 

T x = inf{t > 0, X t > x}. 

Comme X est auto-similaire d'indice 1/a on a T x = x a Ti, et la continuité à droite de X entraîne 
d'autre part {St > x} = {T x < t} p. s. Après changement de variable y = qx a dans Q, on trouve 

(7) / e-^PpTa G dt] = E a (y) - ay 1 ^ E' a (y), y > 0, 

et ceci entraîne que la fonction à droite est complètement monotone avec mesure de Bernstein 
P[Ti G dt]. On va maintenant inverser ([7]) et donner une expression de la loi de T\. Pour a = 2 et 
X = y/2B, le terme de droite vaut et on retrouve l'expression classique pour la densité de T\ 
qui est e _1//4 ' I2t\pwt. On fixe maintenant a G]l,2[ et on note 

fi = inf{t > 0, X t > 1}. 

Il est bien connu et facile à voir - voir le Théorème VII. 1 dans [2\ - que E[e _A ^ 1 ] = e~ Al/a pour tout 
A > 0, autrement dit que T± est une variable stable positive d'indice 1/a. Sa densité peut être 
donnée sous forme intégrale - voir [9] pour le calcul : 

U(t) = - / e~(* u+ " 1/acos ^/ a »sin(u 1 / Q sin(7r/a))d«, 

expression qui se simplifie légèrement pour a = 3/2 en termes de la fonction de Whittaker 
W-l/2,-1/6 " v °i r à nouveau [S]. On considère d'autre part une variable T indépendante de î\ 
et de densité 

= (-^^(l + t 1 /») 
na(t 2 — 2tcos7ra + 1) 

sur M + . Il est évident après un changement de variable dans ([3]) que T est bien une variable aléatoire 
finie p. s. mais je n'ai pas trouvé d'exemples où T jouât un rôle dans la littérature. En particulier, 
malgré les apparences, il ne semble pas que T soit reliée à une loi de l'arcsinus - voir [2j Théorème 
III. 6 - ou encore à un noyau fractionnaire de type Kato - voir jj[2] Formule XI. 11.6, puisque a > 1. 

Théorème 3. Avec les notations précédentes, on a pour tout a G]l,2[ 

Ti = T x Ti. 

Preuve : En combinant (|7|) et le théorème 1, on a pour tout x > 

/ e-^PfTi G dt] = D a (x) = / e~ xt ^i a {dt) 
Jo Jo 
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où l'on a repris la notation ©. D'autre part, comme par définition 



e = l e 

o 



xt I — (xt) a u I a 



f fi (u) du, 



on déduit par Fubini et en faisant le changement de variable s = t a u 



JO JQ 

oo / roo 



ir(t 2a - 2t a cosvra + 1) 
du\ 



ce qui prouve en inversant la transformée de Laplace que T x a pour densité 

Mt) = r f f { U )f T (t/u)- = r Mu)f fi (t/ U )- 

qui est bien celle de T x Ti. □ 

Cette identité en loi peut se décliner de plusieurs façons. On sait par exemple par auto-similarité 
que Ti = S^ a et fi = 5f a , d'où 

(8) 5i = T-^ a x Si 

avec les notations précédentes. Ceci donne une autre relation entre Si et les fonctions de Mittag- 
Leffler que celle donnée par ©. On sait en effet que Si suit une loi de Mittag-Leffler d'indice 1/a 
au sens où 

E[ e - A ^] = Ei /a (-X), A > 

et l'on renvoie pour cela à l'Exercice 29.18 dans [TO] qui donne une preuve probabiliste du résultat 
de Pollard [8]. En utilisant ([8j) et un changement de variable que nous laissons au lecteur, on 
obtient donc une relation probablement inutile, mais inédite à notre connaissance, entre fonctions 
de Mittag-Leffler d'indices différents : 

Corollaire 4. On fixe a g]1, 2[ et on note fi/ a la densité de Bernstein associée à la fonction 
complètement monotone x Eij a (—x). Avec les notations précédentes, on a 

n r \ r f°° , r n -sxVu* (" sin™) (1 + s 1 ^) 

D a (x) = / / /i/ a («)e /u ^-r^ 7T 1 -^duds 

1 1 1 TTa(s z — 2scos7ra + 1) 



J 



pour tout x > 0. 



A l'aide du théorème de Skorohod qui stipule que Si = Xi sachant X\ > - voir pour cela les 
Exercices 29.7 et 29.18 dans [10], on obtient enfin une relation entre les lois de Si et X\\ 

Corollaire 5. Avec les notations précédentes on a 

Si = -(T~ 1/a x Xi) sachant X x < 0. 
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3. UN CALCUL DE CONSTANTE 

On va s'intéresser maintenant au comportement asymptotique de la fonction de répartition de 
T\ en l'infini. On sait par la Proposition VIII. 2 dans 2\ et l'identité T\ = S x a que 

(9) F[T>t] ~ Kt 1 ^' 1 , t^+oo 

où k est une constante dépendant de la loi de la trajectoire du processus X - plus précisément 
de celle de son processus d'échelle bivarié {L~ l ,H) en suivant les notations de [2J. Les résultats 
précédents permettent de calculer k : 

Proposition 6. La constante k dans ([9]) vaut \/T(p)T{\/a). 

Preuve : Dans le cas a = 2 et k = 1/^/tt, le résultat est bien connu et découle directement de 
l'expression de la densité de T\ rappelée plus haut. Dans le cas a 2[, on déduit d'abord de ([7]) 
et du Théorème 1 que pour tout A > 



f Jo \^ u 2a — 2u a cos na + 1 

Après une intégration par parties et un changement de variables, ceci implique 



X a / e - AQ *P[Ti >t]dt = X e- xt ¥[T > t a ] dt 
Jo Jo 



avec les notations précédentes pour T. Grâce à l'expression de fy, on sait d'autre part que 



sin 7ra i , v 

P T > t a ] ; -t 1 ' , t -> OO. 

En utilisant deux fois de suite un théorème taubérien - celui donné par exemple dans [5] Chapitre 
XIII. 5, Théorème 4, on voit que 

prT sin7rar(2-a) l/a _ x l ,, a _ x 

tt{1 — a)T(l/a) r(a)r(l/a) 
où dans la dernière égalité on a utilisé successivement la formule de récurrence et la formule des 
compléments pour la fonction Gamma. □ 

Remarque 7. En écrivant 



A' 



poo roc 

/ e - Aa ¥[Ti <t}dt = X e- M F[T < t a ] dt, A > 0, 
Jo Jo 



on voit à cause de l'asymptotique P[T < t a ] ~ (— sin7ra/7ra;)t Q et à nouveau par un théorème 
taubérien que 

P[Ti<t] ~ -(-r(a)sinvra)t = ( _1 ) t, t^0+, 

résultat déjà connu en combinant l'identité T\ = S x a avec la Proposition VIII. 4 de [2] et l'estimée 
(14.37) de [10] dans le cas particulier p = 1 — 1/a. 
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4. Remarques finales 

En fait le résultat de la Proposition 6 avait lui aussi déjà été obtenu, plus récemment et avec une 
normalisation différente, dans pQ Corollaire 3 Formule (2.61), comme conséquence d'un résultat 
plus général qui était un développement complet en série alternée de la densité de S\ - voir pQ 
Formule (2.8) : 

1 



(10) = 



_ ,cm — l)r(l + lia — n 



-x an ~ 2 . 



Quitte à changer de variable, ce résultat est en effet équivalent à une décomposition du même type 
pour la densité de T\ : 

/TlW = ^aT(an - 1)T(1 + 1/a - n/^ " 

et entraîne donc automatiquement la Proposition 6. Pour obtenir (jlOp . Bernyk, Dalang et Peskir 
avaient inversé la double transformée de Laplace © d'abord en temps et en termes de fonctions 
d'erreur généralisées, puis en espace à l'aide d'une équation intégrale fractionnaire de type Abel, 
dont ils avaient pu donner la solution par un développement en série. Remarquons au passage que 
les fonctions de Mittag-Leffler d'indice plus grand que 1 interviennent dans la résolution d'équations 
intégrales fractionnaires de type Abel [Bj. Mais il ne semble pas que l'on puisse déduire de [6j [T] 
l'identité en loi du Théorème 3. 

La méthode de cette note, où l'on a inversé Q d'abord en espace puis en temps, est plus simple 
que celle de pQ. En revanche, l'expression de la densité de T\ obtenue par le Théorème 3 comme 
une intégrale double : 

, , . -sinTO f°° x t i \ l + u l ' a 

mit) = / f f (t/u)— — —du 

an Jq 1 w 3 — 2u* cos ira + u 

(-sinvra) f 00 / [°° cos(7r/a)) , jl/a f , „ \ 



e - ( tv,u+v> œ SW a)) sin f v l/a ^fa/a)) dv ) , ' " du 

7r z a J \Jq ) u — lu cos ira + u 

ne permet pas de retrouver la décomposition (jlip puisque le terme — v cos(7r/a) > sous 
l'exponentielle interdit d'appliquer Fubini. D'un autre côté, on peut trouver cette dernière ex- 
pression un tantinet plus ramassée que celle obtenue en changeant la variable dans la formule 
(2.56) de pp. 

Une autre formulation de la densité de T\ peut être obtenue en inversant ([7]) comme dans le 
Théorème 1. Avec les mêmes notations on peut en effet écrire 

E[e-^] = E a {x)- ax ^E' a { X ) = ±J j( £^^ \dt. 
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D'autre part, il est possible de choisir H = H$ telle que son intersection avec Re(i) > soit contenue 
dans le cercle centré en l'origine de rayon et son intersection avec Re(i) < dans un secteur 

angulaire — ir/2 — 5 < 6 < ir/2 + 9 avec 5 aussi petit qu'on veut - voir à nouveau [1] Figure 56. 
Comme a G]l,2[ le changement de variable t i— > t a transforme alors H s en un chemin de Hankel 
situé dans le demi-plan Re(t) < x pour ô assez petit. Raisonnant exactement comme dans [9] on 
peut donc calculer 

1 f t a - l é J 1 f e tl/a J 

dt = / dt 

il i - ■> 27rai J H 

oo 



27ri Ju t a — x 2noi\ J u t — x 

« Jo 



-xt 



:t g fi (t)dt 



où l'on a noté 



(12) g (t) = — e^ tu - ul/acos ^hixi{u l / a sm(Tr/a))du. 

Le calcul du deuxième terme se fait comme plus haut et nous laissons les détails au lecteur : 

j_ / ^t dt _ r - u 

27ri J H t a - x i Jo \u 2a - 2u a cos ira + 1 J 

u l / a \ 
f f (t/u)— — — ■ — du dt, 

11 ii** — 2u COS TTOi + U 1 




d'où l'on déduit 

(13) f Tl (t) = (l/a)g fi (t) + (l-l/a)h fi (t) 

en ayant noté 

sm ira /""„,., , / u x 



,1/a 

M*) = 7 ~TT / ffS l l u 



Puisque 



u 2 — 2u cos 7ra + 1 
smna)u l l a 



(a — 1)tï(u 2 — 2ncos7ra + 1) 
est une densité de probabilité sur M + , si l'on note T la v.a. correspondante on voit que hf est la 
densité du produit indépendant T x T\. D'autre part, en comparant (fT3|) avec l'expression de fa 
donnée par le Théorème 3, on voit que est aussi une densité de probabilité, celle du produit 
indépendant T x T\ où T a pour densité 

— sin ira 



ii 



tï{v? — 2ucos7ra + 1) 

En particulier g^, est positive, ce qui n'était pas évident a priori, et la décomposition (fT3"|) se lit 
comme une combinaison convexe de densités de probabilité. Le premier terme (l/a)gf donne le 
comportement de fa en 0, et le second (1 — l/a)hf celui de fa à l'infini. 
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Remarquons enfin que la densité a quelque parenté avec ff - d'où sa souscription en T\ - 
puisque son développement en série alternée, que l'on tire de (fl"2j) exactement comme dans [S], est 



En changeant une dernière fois la variable, on trouve donc l'expression d'une certaine transformée 
intégrale reliant ces deux séries alternées : 



avec laquelle je concluerai, en remerciant Jean-François Burnol et Francis Hirsch pour leur aide 
dans l'élaboration de cet article. 
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tandis que celui de ff est 
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